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The magnetic critical scattering in Sr 2 lr 04 has been characterized using X-ray resonant magnetic 
scattering (XRMS) both below and above the 3D antiferromagnetic ordering temperature, Tjv- The 
order parameter critical exponent below Tn is found to be /3=0.195(4), in the range of the 2D XYhi 
universality class. Over an extended temperature range above Tn, the amplitude and correlation 
length of the intrinsic critical fluctuations are well described by the 2D Heisenberg model with XY 
anisotropy. This contrasts with an earlier study of the critical scattering over a more limited range 
of temperature which found agreement with the theory of the isotropic 2D Heisenberg quantum 
antiferromagnet, developed to describe the critical fluctuations of the conventional Mott insulator 
La 2 Cu 04 and related systems. Our study therefore establishes the importance of XY anisotropy in 
the low-energy effective Hamiltonian of Sr 2 lr 04 , the prototypical spin-orbit Mott insulator. 

PACS numbers: 75.40.-s, 78.70.Ck, 71.70.Ej 


The Ruddlesden-Popper series Sr^+iIr^Oan+i of per- 
ovskite iridates has emerged as a fruitful arena in which 
to explore the effects of electron correlations in the strong 
spin-orbit coupling limit. The first two members of this 
series, single-layer Sr2lr04 (n = 1 ) and bi-layer Sr3lr207 
(n = 2 ), are believed to exemplify a new class of spin- 
orbit Mott insulators. Of central importance to our 
understanding of these materials is the emergence of a 
jefT=l/2 groundstate by the combined action of a strong 
cubic crystal-field and spin orbit interactions on the 
electrons of the Ir^+ ions [H. The weakened electron 
correlations typical of the 5 (i elements then splits the 
jefT=l/2 band, opening a gap, leading to a Mott-like 
state. 

Sr2lr04 in particular has attracted considerable at¬ 
tention because of its striking similarities to La2Cu04 
in terms of both its structural and magnetic properties. 
The magnetic structures and excitations of Sr2lr04 have 
been investigated in a number of X-ray resonant magnetic 
scattering (XRMS) studies [ 11 - 0 ] which have allowed an 
effective low-energy Hamiltonian to be proposed and re¬ 
fined. Sr2lr04 forms an antiferromagnetic structure be¬ 
low Tn 225 K in which the moments are confined to 


the a-b planes and canted to follow rigidly the correlated 
rotation of the oxygen octahedra of the / 4 i/acd crys¬ 
tal structure [^. A resonant inelastic X-ray scattering 
(RIXS) experiment 0 ] has revealed a dispersion relation 
somewhat reminiscent of that displayed by La2Cu04 al¬ 
beit with a lower energy scale and much stronger fur¬ 
ther neighbour couplings, which can be derived from a 
smaller ratio of on-site repulsion over hopping amplitude 
0 . This result suggests that the low-energy isospin dy¬ 
namics of the jefT=l /2 states in Sr2lr04 may, to lead¬ 
ing order, be mapped onto an effective isotropic two- 
dimensional Heisenberg Hamiltonian in agreement with 
predictions by Jackeli and Khaliullin Un¬ 


critical scattering studies provide information com¬ 
plementary to that obtainable from the ordered state. 
For thermally driven transitions in classical systems, is¬ 
sues such as dimensionality, relevant anisotropies, etc., 
can be addressed by determining the critical exponents 
both below and above the transition temperature 11 j. 
One question that naturally arises is: in what ways, if 
any, do the critical fluctuations of a lattice decorated by 
jefT=l/2 isospins differ from the case of S=l /2 spins? For 
La2Cu04 and related 2 D Cu^+ (S=l/ 2 ) systems it has 
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been established that the instantaneous magnetic scatter¬ 
ing function S(Q) is well described by the critical prop¬ 
erties of the 2D quantum S=l/2 Heisenberg antiferro- 
magnet on a square-lattice (2DQHAFSL) [l3 - [l^ . The 
2DQHAFSL model itself has been extensively studied us¬ 
ing a range of theoretical and computational techniques, 
with results that are in broad agreement within the range 
of applicability of the assumptions used Here for 

the sake of brevity we compare our data to the results 
obtained by quantum Monte Carlo (QMC) techniques 


Fujiyama et al. 2l| have measured the critical mag¬ 
netic scattering from Sr 2 lr 04 using XRMS in the inter¬ 
val Tn to Tn -I- 25 K. Their principal results are that 
the critical fluctuations are 2D over the entire temper¬ 
ature interval investigated, with an evolution of the in¬ 
plane correlation length consistent with that expected 
for the 2DQHAFSL model [H The latter was in¬ 
terpreted by Fujiyama et al. as evidence in favour of the 
theoretical proposition by Jackeli and Khaliullin 0 that 
the full isospin Hamiltonian (including anisotropies other 
than that due to Hund’s coupling) can be mapped onto 
an effective, isotropic Heisenberg Hamiltonian. However, 
the good agreement reported between data and predic¬ 
tions of the 2DQHAFSL model required a value of the 
nearest-neighbour exchange coupling of J= 100(10) meV 
which they deduced by fitting their data. A subsequent 
RIXS study 0 of the one-magnon dispersion in Sr 2 lr 04 
provided the value J = 60 meV, and a next-nearest- 
neighbour ferromagnetic coupling J' = — 20 meV, seem¬ 
ingly difficult to reconcile with the data of Fujiyama et 
al. 

Here we report a comprehensive XRMS study of the 
critical magnetic fluctuations in Sr 2 lr 04 for temperatures 
both below and above Tn. Our results provide decisive 
evidence for the importance of XY anisotropy in the ex¬ 
change interactions for Sr 2 lr 04 . Below Tn, accurate val¬ 
ues are obtained for the order parameter /3 revealing that 
Sr 2 lr 04 fits into the 2D XY universality class with ^,4 
anisotropy. Above Tn we greatly extend the tempera¬ 
ture window over which the critical scattering has been 
measured to Tn -I- 73 K, and carefully apply corrections 
for the finite instrumental resolution. We establish that 
in contrast to an earlier experimental study for Sr 2 lr 04 
21j the intrinsic isospin fluctuations for T > Tn -|- 5 K 
are dominated by XY type anisotropy. 

The critical scattering experiments were performed on 
beamline 116, Diamond Light Source, UK, and P09, PE¬ 
TRA HI, DESY, Germany. These experiments exploited 
the large enhancement of the X-ray resonant magnetic 
scattering cross-section at the Ir L 3 edge. The sin¬ 
gle crystals of Sr 2 lr 04 used in this study (dimensions 
2 X 1X 0.05 mm^) were flux grown from phase-pure poly¬ 
crystalline Sr 2 lr 04 using techniques described elsewhere 
22l |. and attached to the copper sample mount of a 
closed-cycle refrigerator. This was in turn mounted on a 
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FIG. 1. Order parameters of the magnetic critical scattering 
in Sr 2 lr 04 : (a) integrated intensity and power law fit (solid 
line) corresponding to best-fit values of rN=226.80(6) K and 
/I=0. 195(4). (b): Goodness of fit (purple squares) and 

resulting exponent P (green circles) as function of chosen or¬ 
dering temperature Tm- The fit in (a) corresponds to the 
minimum in x^- 


six-circle diffractometer configured to operate in a verti¬ 
cal scattering geometry. The energy of the incident pho¬ 
ton beam was set to 11.218 keV, just below the T 3 edge of 
iridium, a value found to maximise the intensity of the X- 
ray resonant magnetic scattering. The incident beam size 
was determined to be 200x20 iirn? (HxV). The polariza¬ 
tion of the scattered X-rays was determined by using a Au 
(333) crystal analyser mounted on the detector arm. The 
temperature was measured to a precision of ±0.01 K via 
a thermocouple secured to the sample mount by Teflon 
tape. The wavevector resolution of the instrument, in¬ 
cluding the effects of sample mosaic, was determined by 
mapping Bragg peaks in reciprocal space and was found 
to be typically 1 . 1 x 10 “^ and 3.7xl0“^A“^ perpendicu¬ 
lar and parallel to Q in the scattering plane respectively, 
and 1.4 x 10“^A“^ out of the plane. 

The first objective was to determine Tn and /3 from the 
temperature dependence of the magnetic peak intensity 
Im below Tn. The results are summarized in Fig. [T] A 
fit to the simple power-law form Im oc [(Tn — T)/Tn]^^ is 
shown in Fig. [TJa) which yields /3 = 0.195(4). This value 
is consistent with that provided by neutron scattering 
measurements [ 2 ^, [ 2 ^ . The value of /3 deduced by this 
analysis deviates signihcantly from theoretical values for 
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FIG. 2. Critical scattering data for T > Tn obtained in scans 
performed perpendicular to the wavevector transfer, Q, of the 
appropriate antiferromagnetic Bragg peak. Solid lines: fits 
of a Lorentzian line shape convoluted with the experimental 
resolution. The horizontal bars nnder the peaks represent the 
resolution full width at half maximum. 


both 2D Ising (/3 = 1/8) and 3D systems (/3 ^ 0.35), but 
rather is consistent with value for the 2D XYh^ univer¬ 
sality class in the strong anisotropy limit [25| . 

Next, a detailed investigation of the critical scattering 
above Tn was undertaken. For these studies the high 
photon energy and relatively broad energy resolution of 
X-ray diffractometers (> 1 eV) offer an advantage over 
neutrons in that they provide an accurate frequency in¬ 
tegration to yield the instantaneous magnetic scattering 
function S(Q). On the other hand the intrinsic high 
wavevector resolution of X-ray techniques presents a chal¬ 
lenge in terms of following weak critical magnetic scat¬ 
tering to high temperatures as it broadens and weakens 
further. X-ray experiments to determine the magnetic 
critical scattering above Tn have revealed that the crit¬ 
ical fluctuations just above Tn have two components: a 
sharp “central” peak, typically with a Lorentzian squared 
lineshape, believed to be an extrinsic feature due to the 
presence of defects; and a broader, weaker peak with a 
Lorentzian line shape arising from intrinsic critical fluc¬ 


tuations [llj, l2g . Realisation that the magnetic critical 
scattering above Tn can have two such components mir¬ 
rors earlier results for structural phase transitions 27- 


Representative scans of the magnetic critical scatter¬ 
ing above Tn are shown in Fig. [2l The critical scattering 
could be followed out to Tn -I- 73 K, roughly trebling the 
region probed above Tn in Ref. [2]|, which, as we shall 
show, places a much tighter constraint on the relevant 
form of the effective Hamiltonian. The critical scattering 



FIG. 3. Temperature dependence of the correlation length 
(top panel) and the amplitude (bottom panel) of the crit¬ 
ical scattering for Sr 2 lr 04 as filled circles. Open squares 
are data taken from Fujiyama [2^ divided by no. QMG re¬ 
sults with J = 60 meV (100 meV) as a dashed line (dot- 
dashed line) Ref. [T^. Solid line: fit of 2D model with 
XY anisotropy as discussed in text. Inset: Linear scaling 
of ln(^5'o) ~ (3 — 77 ) ht~^. 


above Tn was fitted with several different lineshapes con¬ 
voluted with the instrumental resolution function. The 
most satisfactory over most of the temperature range was 
a Lorentzian, as expected for intrinsic magnetic critical 
scattering. However for Tn < T < Tn -|- 5.5 K a sharper, 
second component appeared to develop, presumably due 
to extrinsic, defect mediated scattering. As it proved dif¬ 
ficult to obtain unique fits in this temperature interval 
we do not further consider the data taken close to Tn. 


In Fig. [3 we compare our data with with various the¬ 
ories and earlier experimental results [^ . The broken 
lines represent QMC results for the 2DQHAFSL model 
for which ^(T) = 0.276 and Sq oc ^^T^, 

where uq = 3.9 A is the Ir-Ir nearest neighbour dis¬ 
tance. It is evident that in the 2D regime the values 
of ^ and So of the intrinsic Lorentzian component have 
much stronger temperature dependences than the expec¬ 
tations of the 2DQHAFSL. The 2DQHAFSL model can¬ 
not therefore describe the critical fluctuations in Sr 2 lr 04 
whatever value of J is chosen. This important realisation 
disagrees with the main conclusion of Ref. 211. 

We therefore considered alternative models. In the 2D 
XY model a phase transition occurs at the Kosterlitz- 
Thouless temperature, Tkt [3I|- Above Tkt the correla¬ 
tion length ^ and scattering amplitude So are described 
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by: ^ oc exp and Sq oc exp [b {2 — rf )/, where b 

is a non-universal constant (typically b^i 1.9), 77 = 1/4 and 
t= (T — Tkt) / Tkt the reduced temperature [S^- For 
real systems, three-dimensional order at a temperature 
2 n > 2 kt occurs as a result of weak interlayer coupling 
J' amplified by the diverging size of correlated regions 
^(T)^. However attempts to fit the pure 2D XY model 
to our data yielded unphysical values of b = 3.0(7) and 
J'/J = 5x 10“®, and were thus discarded. 

Instead, a consistent description of our data was 
obtained using the 2D anisotropic Heisenberg model 
(2DAH), where leading Heisenberg interactions are aug¬ 
mented by easy-plane, XY anisotropy [s^. The 2DAH 
model also has a correlation length that diverges expo¬ 
nentially towards a finite temperature Tkt, placing it in 
the 2D XY universality class. In the limit of long cor¬ 
relation lengths, ^ and Sq scale the same way as for the 
pure 2D XY model. However with increasing tempera¬ 
ture there is a crossover towards ^(T) = ,^0 exp [b/t) and 
Sq oc exp [6 (2 — 77 )//] for ^ < 100 oq [s^. Our data fall 
in this regime. 

The solid lines in Fig.[3]show the fit to the 2DAH which 
is seen to provide an excellent description of the data in 
the 2 D regime. From the fits we obtain ^0 = 0.9(1) A, 
Tkt = 162(10) K and 6 = 2.1(7), with the latter in good 
agreement with theory and experiment for other layered 
materials [sll, [s^, . From the value of the 2D correla¬ 

tion length at Tm it is possible to estimate the ratio of 
the in-plane and inter-layer couplings J'/ J ~ = 

exp(—26/tAr) K. 0.001. This ratio is an order of magni¬ 
tude higher than for the cuprates [mu , which can be 
understood by comparing the shapes of the ground state 
orbitals of Sr 2 lr 04 and La 2 Cu 04 . The jeff=l/2 state is 
approximately cubic, whereas for La 2 Cu 04 the dx 2 _y 2 
orbital lies predominantly in the xy plane. Thus larger 
out-of-plane interactions would be expected for Sr 2 lr 04 
compared with La 2 Cu 04 . The description of ^{T) by in¬ 
clusion oi XY anisotropy is consistent with the extracted 
value of /3 = 0.195(4), which places the fluctuations be¬ 
low T^v in the 2D AT/ 7.4 universality window, where a 
4-fold in-plane anisotropy pushes 0 from 0.231 for 2D 
AT towards 0 = 1/8 for 2D Ising [2^ interactions. 

In general, the effect of AT anisotropy on the spin 
wave dispersion is to lift the degeneracy between the 
Heisenberg modes. The in-plane fluctuations a;|| main¬ 
tain a continuous symmetry and are hence gapless (Gold- 
stone mode), whereas the out-of-plane fluctuations uj± 
are gapped. At the zone centre the out-of-plane spin 
wave gap is given by loq = 4 JS'Tc-\/2Aa, where 0 < Aa < 1 
parametrizes the easy-plane anisotropy 36| . In data from 
a recent high-resolution RIXS experiment on Sr 2 lr 04 
there appears to be a magnon gap of ~30 meV at 
the magnetic zone centre, qualitatively consistent with 
the existence of AT anisotropy. The degree of AT 
anisotropy required to produce a gap of this size ma’ 


be determined by extending the model of Kim et al. p 


to include anisotropic exchange. The resulting 2DAH 
Hamiltonian reads 

nab = J2J [S/SJ + SfS^ + (1 - Aa) S/S^] 

0,j) 

+ ^ J2S^Sj + •^3S^SJ, ( 1 ) 

(dj)> ((bj»> 

Here J = Jiso/ (1 — Aa) is the effective nearest-neighbour 
(nn) exchange parameter, Jj^o is the isotropic Heisenberg 
nn exchange, and J 2 , J 3 symbolizes the exchange between 
next-nearest and third nearest neighbours respectively. 
Out-of-plane exchange coupling has been neglected since 
as demonstrated above, it is a factor of ^ 10^ weaker than 
the in-plane terms, and thus not resolvable with RIXS at 
present. 
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FIG. 4. Fit of the anisotropic 2D Heisenberg model to the 
experimental single magnon dispersion. Squares: Energy (a) 
and spectral weight (b) of the single magnon peak obtained 
from fitting RIXS data provided in Ref. [^. Errorbars are 
smaller than the displayed symbols. Solid and dashed lines: 
best fit of gapped and gapless modes to experimental data 
with J = 57(7) meV, J 2 = —18(3) meV, J 3 = 14(2) meV, 
Aa = 0.08(1). 

We fitted the theoretical dispersion of this 2DAH 
model to the experimental data from Ref. @ using a 
linear spin-wave model derived from Eq. [1] with the re¬ 
sults shown in Fig. 21 There is good agreement between 
the experimental dispersion and that calculated from the 
anisotropic model, with the fitted exchange constants es¬ 
sentially identical to those obtained from the earlier low- 
resolution study within the error bars [H . The agreement 
for the spectral weight is perhaps less compelling but 
can be said to be qualitatively reasonable, in particular 
the model reproduces the non-zero spectral weight at the 
zone centre. Moreover, it proved difficult to extract the 
spectral weight from the RIXS data as there appears 
to be a significant variation of lineshape as a function of 
momentum transfer, which may be evidence for an ad¬ 
ditional excitation mode between the single magnon and 
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bi-magnon peaks (the latter was accounted for in the fit¬ 
ting procedure). Fitting the dispersion and intensities to 
the localized spin model proposed by Igarashi [37| pro¬ 
vided similar results, corroborating their prediction of an 
out-of-plane gap for Sr 2 lr 04 . The easy-plane anisotropy 
parameter Aa = 0.08(1) determined here should be con¬ 
sidered as an upper bound to the true anisotropy, since 
the finite momentum and energy resolution means that 
the excitation energy at the zone centre cannot be fully 
resolved. Nevertheless this is significantly larger than 
that observed for La 2 Cu 04 (Aa = 2.0(5) x 10“^) [s^ . 
which further illustrates the relative importance oi XY 
anisotropy for Sr 2 lr 04 . 


In conclusion, we have presented a detailed study of 
the critical magnetic properties of the single-layer per- 
ovskite iridate Sr 2 lr 04 . Our results establish the im¬ 
portance of XY anisotropy in its effective low-energy 
Hamiltonian. Although in their seminal paper Jack- 
eli and Khaliullin [l^ proposed that the interactions in 
Sr 2 lr 04 are isotropic (Heisenberg) they note that the 
leading anisotropy term should have XY character. Here 
we have shown that the XY anisotropy is sufficiently 
strong to push the critical properties of Sr 2 lr 04 signifi¬ 
cantly away from the isotropic 2DQHAFSL model. Thus 
while Sr 2 lr 04 and La 2 Cu 04 share some similarities, as 
far as their effective magnetic interactions and critical 
properties are concerned, they are quite distinct. 
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